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Rapid  Computation  of  Yield  Tables  for 
Managed,  Even-Aged  Timber  Stands 

Clifford  A.  Myers  and  Gary  L.  Godsey 


Introduction 

Yield  tables  for  managed  stands  are  valuable 
guides  for  those  who  determine  and  for  those 
who  carry  out  timber  management  policies. 
Reported  volumes  are  probable  yields  that  result 
from  specified  combinations  of  such  factors  as 
frequency  and  intensity  of  thinning,  site  quality, 
and  utilization  standards.  The  volumes  may  be 
used  to  compute  various  money  yields  and  costs 
(Myers  1968).  There  is,  however,  the  disad- 
vantage that  a  yield  table  reports  probable 
yields  for  only  one  combination  of  stand  char- 
acteristics and  management  decisions.  Managers 
interested  in  comparing  outcomes  of  several 
of  the  many  alternatives  possible  do  not  find 
a  yield  table  to  be  an  effective  basis  for 
decisions.  But  a  group  of  yield  tables  can  be 
a  useful  tool  for  decisionmaking  so  long  as  the 
tables  are  not  too  costly  in  time  and  money. 

Groups  of  yield  tables  can  be  obtained  from 
measurements  of  trees  on  temporary  plots  and 
computations  made  on  digital  computers.  Use 
of  temporary  instead  of  permanent  plots  avoids 
(1)  long  delays  before  yield  tables  become 
available,  (2)  possible  loss  of  expensive  in- 
stallations, and  (3)  the  possibility  that  changes 
in  utilization  and  management  practices  will 
render  the  plots  obsolete  before  much  infor- 
mation can  be  obtained.  Comparisons  of  data 
from  temporary  and  permanent  plots  show  that 
data  from  temporary  plots  are  entirely  adequate, 
and  permit  use  of  several  approaches  to  yield 
table  construction  (Decourt  1965,  Myers  1966, 
Vuokila  1965).  Once  several  relationships  be- 
tween stand  variables  have  been  established, 
a  manager  can  examine  the  probable  outcomes 
of  many  possible  variations  in  management. 
There  is  no  need  to  delay  decisions  or  to 
speculate  on  what  may  happen  if  a  condition 
or  procedure  changes.  Large  numbers  of  tables, 
each  based  on  a  specified  set  of  alternatives, 
can  be  computed  and  printed  at  a  cost  of  a  few 
cents  each. 


Procedure  for  computation  of  yield  tables  is 
presented  here  as  two  computer  programs 
written  in  FORTRAN  IV.  Measurements  needed 
and  their  use  in  mathematical  expression  of 
useful  relationships  are  described  in  detail  else- 
where (Myers  1966,  1967).  Parts  of  this  material 
are  repeated  in  the  section  headed  "Information 
Needed."  The  first  program  described, program 
YIELD  (appendix  1),  computes  and  prints  a 
yield  table  for  each  specified  combination  of 
alternatives.  Program  THIN1  (appendix  2) 
calculates  the  increases  in  stand  diameter 
caused  by  various  intensities  of  initial  thinning 
from  below.  These  data  are  used  to  calculate 
a  relationship  used  in  program  YIELD.  Defini- 
tions of  variables  used  in  each  program  appear 
at  the  beginning  of  the  program. 

The  two  programs  and  test  problems  (ap- 
pendices 1  and  2)  use  as  examples  relationships 
established  for  ponderosapine  (Pinuiponderosa 
Laws.)  in  the  Black  Hills  of  South  Dakota  and 
Wyoming.  Both  programs  may  be  converted 
easily  to  apply  to  other  standards  of  utilization 
and  management,  or  to  other  tree  species. 


Information  Needed 


Computation  of  yield  tables  from  temporary 
plot  data  requires  at  least  seven  items  of  in- 
formation. Additional  items  may  be  necessary 
for  some  species.  For  example,  mortality  may 
be  of  importance  in  thinned  stands,  although  it 
was  not  for  the  ponderosa  pine  used  as  the 
example.  The  seven  basic  items  needed  and 
procedures  to  obtain  them  are: 

1.  Description  of  unthinned,  young  stands. 

Stands  with  the  ages  and  densities  expected 
at  possible  times  of  initial  thinning  are  sampled. 
Age  at  initial  thinning,  number  of  trees  per 
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acre,  and  average  stand  diameter  are  obtained 
for  use  as  initial  entries  in  yield  tables.  A 
variety  of  stands  is  sampled  if  age  and  density 
at  time  of  initial  thinning  are  variables  of 
interest.  Each  plot  also  provides  a  cumulative 
stand  table  with  numbers  of  trees  in  percent 
(Bruce  and  Schumacher  1950),  for  use  in  com- 
putation of  postthinning  diameters  by  program 
THIN1.  Regressions  of  average  diameter  on 
density  and  age,  and  of  percentage  of  trees  on 
average  diameter  are  established  graphically  or 
mathematically  to  obtain  data  on  intermediate 
stand  conditions  not  sampled. 


2.  Diameter  increase  from  growth. 

Sample  plots  in  stands  that  represent  a 
wide  range  of  densities,  average  diameters,  and 
site  qualities  are  measured  to  obtain  a  means 
of  estimating  periodic  diameter  growth.  Plot 
trees  are  bored  at  breast  height  to  determine 
radial  wood  growth.  Data  are  obtained  to  relate 
diameter  inside  bark  to  diameter  outside  bark 
and  vice  versa.  Periodic  wood  and  bark 
growths  are  subtracted  from  present  diameters 
to  get  diameters  at  the  beginning  of  the  period. 
Present  and  past  diameters  are  grouped  in 
stand  tables,  with  diameters  of  recently  dead 
trees  added  to  the  past  stand  tables.  Present 
and  past  stand  tables  are  used  to  compute 
present  and  past  basal  areas  per  acre  and  aver- 
age stand  diameters. 

Appropriate  values  from  the  past  stands  are 
used  as  independent  variables  to  estimate  aver- 
age diameters  at  the  end  of  the  projection 
period.  For  example,  future  diameters  of  pon- 
derosa  pine  can  be  estimated  from  initial 
diameter,  initial  basal  area,  and  site  index. 
The  equation  for  a  10-year  projection  period 
appears  in  program  YIELD  as  the  long  FOR- 
TRAN statement  for  DBHO. 

The  sample  plots  can  also  supply  other 
needed  values,  such  as  heights  and  a  stand 
volume  equation. 


3.  Diameter  increase  from  thinning. 

Change  in  average  stand  diameter  caused 
by  intermediate  cuttings  subsequent  to  initial 
shinning  can  be  estimated  by  trial  marking  of 
temporary  plots.  An  average  increase  of  0.4 
inch  is  used  in  the  example  of  program  YIELD. 
Rethinning  to  about  the  same  density  level 


(item  4)  each  time  is  assumed.  Amount  of 
increase  was  not  correlated  with  tree  size  or 
stand  density  for  the  stands  examined. 

A  method  of  estimating  the  increase  in 
diameter  caused  by  initial  thinning  is  given  in 
program  THIN1. 


4.  Stocking  guide. 

A  guide  to  the  growing  stock  to  be  left 
after  each  intermediate  cutting  is  prepared  as 
a  relationship  between  basal  area  and  average 
stand  diameter  (table  1).  Results  of  thinning 
studies  and  data  from  temporary  plots  are  used 
to  establish  a  graph  of  desirable  basal  area  over 
diameter  for  local  average  site  quality.  "Best" 
stand  density  for  each  average  diameter  is  based 
on  such  criteria  as  production  in  cubic  feet  and 
probable  length  of  saw  log  rotations.  Desired 
stand  density  varies  with  the  objectives  of 
management,  so  the  series  of  "best"  densities 
must  be  converted  to  a  general  relationship. 

The  curve  of  basal  area  on  diameter  can  be 
considered  as  a  guide  curve  from  which  a 
family  of  curves  can  be  drawn  or  computed. 
Basal  areas  of  table  1  increase  with  diameter 
until  10.0  inches  is  reached,  and  remain  constant 
thereafter.  Basal  area  is  80  square  feet  at  10.0 
inches  diameter,  and  the  series  is  labeled  level 
80.  Other  possible  series  are  identified  by  the 
basal  area  desired  when  average  stand  diameter 
is  10.0  inches  or  greater.  Basal  areas  of  any 
stocking  level  can  be  computed  by  multiplying 
the  values  for  level  80  (table  1)  by  the  ratio 
level/80.  For  example,  basal  areas  for  level 
100  are  each  100/80  times  the  corresponding 
basal  areas  of  table  1. 

Any  desired  form  of  the  guide  curve  for 
the  basal  area-diameter  relationship  may  be 
used.  It  will  be  necessary  to  modify  appropriate 
statements  of  program  YIELD  for  other  bases 
than  level  80. 


5.  Dominant  and  codominant  heights. 

Measurements  on  sample  plots  or  values 
from  good  site  index  curves  provide  height-age- 
site  index  relationships  (table  2).  Average 
heights  of  dominant  and  codominant  trees  are 
usually  more  useful  than  other  average  heights. 
They  are  changed  less  than  other  averages 
when  stands  are  cut  according  to  most  guides 
for  intermediate  cutting. 
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Table  1. — Basal  areas  after  intermediate  cutting  in  relation  to  average  stand  diameter 
Growing  stock  level  80  for  Black  Hills  ponderosa  pine 


Average 

Basal 

Average 

Basal 

Average 

Basal 

Average 

Basal 

stand  d.b.h. 

area 

stand  d.b.h. 

area 

stand  d.b.h. 

area 

stand  d.b.h. 

area 

after  cutting 

per 

after  cutting 

per 

after  cutting 

per 

after  cutting 

per 

(Inches) 

acre 

( Inches ) 

acre 

(Inches) 

acre 

(Inches) 

acre 

Sq.  ft. 

Sq.  ft. 

Sq.  ft. 

Sq.  ft. 

2.0 

12.3 

4.0 

35.2 

6.0 

56.6 

8.0 

72.5 

2.1 

13.3 

4.1 

36.4 

6.1 

57.6 

8.1 

73.2 

2.2 

14.4 

4.2 

37.6 

6.2 

58.5 

8.2 

73.8 

Z.  .  J 

15.5 

4.3 

38 . 7 

6 . 3 

D  7  ■  H 

O  .  J 

74.4 

2.4 

16.6 

4.4 

39.8 

6.4 

60.3 

8.4 

74.9 

2.5 

17.7 

4.5 

41.0 

6.5 

61.1 

8.5 

75.4 

2.6 

18.9 

4.6 

42.2 

6.6 

62.0 

8.6 

75.8 

2.7 

20.1 

4.7 

43.4 

6.7 

62.9 

8.7 

76.3 

2 . 8 

21 . 3 

4 . 8 

44 . 6 

6 . 8 

63 . 8 

8 . 8 

7fi  7 

2.9 

22.5 

4.9 

45.7 

6.9 

64.7 

8.9 

77.1 

3.0 

23.6 

5.0 

46.7 

7.0 

65.5 

9.0 

77.5 

3.1 

24.8 

5.1 

47.7 

7.1 

66.2 

9.1 

77.9 

3.2 

26.0 

5.2 

48.8 

7.2 

67.0 

9.2 

78.2 

3.3 

27.2 

5.3 

49.9 

7.3 

67.7 

9.3 

78.5 

3.4 

28.4 

5.4 

50.9 

7.4 

68.4 

9.4 

78.8 

3.5 

29.6 

5.5 

51.8 

7.5 

69.1 

9.5 

79.0 

3.6 

30.7 

5.6 

52.8 

7.6 

69.8 

9.6 

79.2 

3.7 

31.8 

5.7 

53.7 

7.7 

70.5 

9.7 

79.5 

3.8 

32.9 

5.8 

54.7 

7.8 

71.2 

9.8 

79.7 

3.9 

34.1 

5.9 

55.6 

7.9 

71.9 

9.9 

79.9 

10.0+  80.0 


Table  2. — Average  height  of  dominant  and  codominant  trees  at 
various  ages,  Black  Hills  ponderosa  pine 


Main 
stand  age 

Site  index 

class 

(Years) 

40 

1  50 

|  60 

70 

-  -  -  -  -  Feet 

10 

4.5 

4.5 

4.5 

4.5 

20 

9 

10 

12 

16 

30 

11 

16 

20 

26 

40 

17 

22 

28 

35 

50 

21 

28 

35 

43 

60 

26 

33 

41 

50 

70 

30 

38 

47 

56 

80 

34 

43 

52 

61 

90 

37 

47 

57 

66 

100 

40 

50 

60 

70 

110 

43 

53 

63 

74 

120 

45 

56 

66 

77 

130 

46 

59 

69 

80 

140 

48 

61 

71 

83 

150 

50 

63 

73 

86 

160 

51 

64 

75 

88 

170 

52 

65 

77 

90 

180 

53 

66 

78 

91 
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6.  Stand  volume  equation. 


7.  Volume  conversion  factors. 


Volumes  measured  on  temporary  growth 
plots  and  elsewhere  provide  the  data  needed 
to  compute  stand  volume  equations.  Multiple 
regression  methods  are  used  to  compute  an 
equation  for  predicting  total  cubic  feet  per 
acre  from  other  common  stand  measures.  In- 
dependent variables  most  useful  for  ponderosa 
pine  stands  are  basal  area,  average  height  of 
dominant  and  codominant  trees,  and  average 
stand  diameter.  The  equation  app  ears  at  various 
places  in  program  YIELD  as  statements  for 
TOTO  and  TOTT. 

Table  3. — Factors  for  conversion  of  stand  volumes  in  total  cubic  feet  to 
merchantable  cubic  feet  per  acre,     Black  Hills  ponderosa  pine 


Stand  volumes  in  total  cubic  feet  can  be 
converted  to  volumes  in  cubic  or  board  feet 
for  any  standard  of  utilization.  Plot  volumes 
in  each  measure  and  standard  of  interest  are 
compared  with  corresponding  total  cubic 
volumes  to  obtain  ratios  or  conversion  factors 
(tables  3  and  4).  Volume  standards  may  be 
changed  without  change  in  the  basic  set  of 
computations  that  produce  yield  tables. 


Average 

Average 

Average 

stand 

Ratio 

of 

stand 

Ratio 

of 

stand 

Ratio  of 

diameter 

merchantab le 

diameter 

merchantable 

diameter 

merchantable 

(Inches) 

to  total 

volume 

(Inches) 

to  total 

volume 

(Inches) 

to  total  volume 

5.0 

0. 332 

8.1 

0.849 

11. 

9 

0.940 

5.1 

.355 

8.2 

.856 

12. 

1 

.941 

5.2 

.377 

8.3 

.862 

12. 

4 

.942 

5.3 

.400 

8.4 

.868 

12. 

7 

.943 

5.4 

.422 

8.5 

.872 

12. 

9 

.944 

5.5 

.444 

8.6 

.876 

13. 

1 

.945 

5.6 

.465 

8.7 

.880 

13. 

3 

.946 

5.7 

.487 

8.8 

.884 

13. 

5 

.947 

5.8 

.508 

8.9 

.888 

13. 

7 

.948 

5.9 

.530 

9.0 

.892 

13. 

9 

.949 

6.0 

.552 

9.1 

.896 

14. 

2 

.950 

6.1 

.575 

9.2 

.899 

14.4 

.951 

6.2 

.597 

9.3 

.902 

14. 

7 

.952 

6.3 

.618 

9.4 

.906 

14. 

9 

.953 

6.4 

.639 

9.5 

.910 

15. 

2 

.954 

6.5 

.659 

9.6 

.913 

15.4 

.955 

6.6 

.678 

9.7 

.916 

15. 

8 

.956 

6.7 

.694 

9.8 

.920 

16. 

3 

.957 

6.8 

.710 

9.9 

.923 

16. 

8 

.958 

6.9 

.725 

10.0 

.926 

17. 

3 

.959 

7.0 

.740 

10.1 

.928 

17. 

8 

.960 

7.1 

.753 

10.2 

.930 

18. 

3 

.961 

7.2 

.766 

10.3 

.931 

18. 

8 

.962 

7.3 

.778 

10.4 

.932 

19. 

3 

.963 

7.4 

.789 

10.5 

.933 

19. 

8 

.964 

7.5 
7.6 
7.7 
7.8 
7.9 
8.0 


,799 
,809 
,818 
,826 
.834 
.842 


10.7 
10.9 
11. 1 
11.3 
11.5 
11.7 


,934 
.935 
,936 
.937 
.938 
.939 


20.3 
20.9 
21.7 
22.5 
23.3 
23.9 


965 
966 
,967 
968 
,969 
,969 


To  4.0- inch  top  in  trees  6.0  inches  d.b.h.  and  larger. 

Factor  for  an  unlisted  diameter  equals  factor  for  next  smaller  listed 
diameter.     For  example,  factor  for  15.6  inches  is  .955. 
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Table  4. — Factors  for  conversion  of  stand  volumes  in  total  cubic  feet  to  board  feet 
Scribner  rule  per  acre,1  Black  Hills  ponderosa  pine 


Ratio  of 

Average 

Ratio  of 

Average 

Ratio  of 

Average 

Ratio  of 

stand 

board  feet 

stand 

board  feet 

stand 

board  feet 

stand 

board  feet 

diameter 

to 

diameter 

to 

diameter 

to 

diameter 

to 

cubic  feet 

( Inches ) 

cubic  feet 

(Inches) 

cubic  feet 

(Inches ) 

cubic  feet 

8.0 

0.78 

11.9 

3.49 

15.8 

4.71 

19.7 

5.42 

8.1 

.85 

12.0 

3.56 

15.9 

4.73 

19.8 

5.44 

8.2 

.92 

12.1 

3.61 

16.0 

4.76 

19.9 

5.45 

8.3 

.99 

12.2 

3.65 

16.1 

4.78 

20.0 

5.46 

8.4 

1.06 

12.3 

3.69 

16.2 

4.81 

20. 1 

5.47 

8.5 

1.13 

12.4 

3.  73 

16.3 

4.83 

20.  2 

5 . 48 

8.6 

1.20 

12.5 

3.77 

16.4 

4.86 

20.  3 

5 . 50 

8.7 

1.27 

12.6 

3.80 

16.5 

4.88 

20.4 

5.51 

8.8 

1.34 

12.7 

3.84 

16.6 

4.90 

20.5 

5.52 

8.9 

1.41 

12.8 

3.88 

16.7 

4.92 

20.  6 

5.53 

9.0 

1.48 

12.9 

3.91 

16.8 

4.94 

20.7 

5.54 

9.1 

1.55 

13.0 

3.95 

16.9 

4.96 

20.8 

5.56 

9.2 

1.62 

13.1 

3.98 

17.0 

4.98 

20.9 

5.57 

9.3 

1.68 

13.2 

4.02 

17.1 

5.00 

21.0 

5.58 

9.4 

1.75 

13.3 

4.05 

17.2 

5.02 

21. 1 

5.59 

9.5 

1.82 

13.4 

4.08 

17.3 

5.04 

21.  2 

5 . 60 

9.6 

1.89 

13.5 

4.11 

17.4 

5.06 

21.3 

5.  61 

9.7 

1.96 

13.6 

4.14 

17.5 

5.08 

21.4 

5.62 

9.8 

2.03 

13.7 

4.17 

17.6 

5.10 

21.5 

5.63 

9.9 

2.10 

13.8 

4.20 

17.7 

5.12 

21.  6 

5 . 64 

10.0 

2.17 

13.9 

4.23 

17  .8 

5.13 

21.  7 

5 . 65 

10.1 

2.24 

14.0 

4.25 

17.9 

5.15 

21.8 

5 .  66 

10.2 

2.31 

14.1 

4.28 

18.0 

5.17 

21.9 

5.67 

10.3 

2.38 

14.2 

4.31 

18.1 

5.19 

22.0 

5.68 

10.4 

2.45 

14.3 

4.34 

18.  2 

5.  21 

22.1 

5.  69 

10.5 

2.52 

14.4 

4.37 

18.3 

5.22 

22.  2 

5.70 

10.6 

2.59 

14.5 

4.39 

18.4 

5.  24 

22.3 

5.71 

10.7 

2.65 

14.  6 

4.42 

18.5 

5.26 

22.4 

5.72 

10.8 

2.72 

14.7 

4.44 

18.6 

5.27 

22.5 

5.73 

10.9 

2.79 

14.8 

4.47 

18.7 

5.29 

22.6 

5 . 74 

11.0 

2.86 

14.9 

4.49 

18.8 

5.30 

22.7 

5.75 

11.1 

2.93 

15.0 

4.52 

18.9 

5.32 

22.8 

5.76 

11.2 

3.00 

15.1 

4.54 

19.0 

5.33 

22.9 

5.77 

11.3 

3.07 

15.2 

4.56 

19.1 

5.35 

23.0 

5.78 

11.4 

3.14 

15.3 

4.58 

19.2 

5.36 

23.1 

5.79 

11.5 

3.21 

15.4 

4.61 

19.3 

5.37 

23.2 

5.80 

11.6 

3.28 

15.5 

4.64 

19.4 

5.39 

23.3 

5.81 

11.7 

3.35 

15.6 

4.66 

19.5 

5.40 

23.4 

5.82 

11.8 

3.42 

15.7 

4.68 

19.6 

5.41 

23.5 

5.83 

To  8-inch  top  in  trees  10.0  inches  d.b.h.  and  larger. 
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Description  of  Program  YIELD 

Operations  performed  by  program  YIELD 
are  identified  by  the  comment  statements  of 
the  source  program  (appendix  1).  Initial  stand 
conditions  and  other  values  are  read  into  com- 
puter memory  in  the  order  and  format  given 
in  the  tabulation  of  order  and  contents  of  the 
data  deck.  The  number  of  yield  tables  com- 
puted and  printed  is  determined  by  the  values 
assigned  NTSTS  on  card  type  1  and  MIX  on 
card  type  2.  NTSTS  is  the  number  of  sets  of 
tables  produced.  MIX  is  the  number  of  tables 
in  a  set  or  the  number  of  stocking  levels 
(DSTY)  tested.  Growing  stock  level  for  initial 
thinning  may  differ  from  that  specified  for 
rethinnings. 

Diameter  after  initial  thinning  is  estimated 
from  diameter  before  thinning  and  an  estimate 
of  the  percentage  of  trees  to  be  retained.  Com- 
putations are  listed  in  FORTRAN  statements 
8  to  13.  The  first  estimate  of  diameter,  DBHE, 
is  tested  against  the  values  in  or  computed 
from  table  1  for  the  specified  growing  stock 
level.  The  estimate  is  modified  as  necessary 
until  diameter  and  basal  area  of  the  yield  table 
and  of  the  growing  stock  level  are  equal.  De- 
rivation of  the  equation  for  PDBHE  (logarithm 
of  DBHE)  is  explained  in  the  description  of 
program  THIN1. 

Clearcutting  at  the  oldest  age  of  interest  is 
assumed.  The  program  can  be  modified  to 
show  the  reservation  of  volume  in  a  shelter- 
wood  or  seed  trees,  if  desired.  Present  as- 
sumptions apply  to  many  species.  For  others, 
a  change  may  be  an  unnecessary  refinement 
where  length  of  the  period  necessary  for  natural 
regeneration  has  not  yet  been  established. 

Replacement  of  several  statements  and  data 
cards  will  modify  the  program  for  other  species. 
Replacements  needed  are: 

1.  Tables  2  to  4,  inclusive.  It  may  be  necessary 
to  change  the  DIMENSION  statement  and 
statements  that  compute  subscripts  for  table 
look-up. 

2.  Table  1,  if  desired.  New  statements  for 
BREAK  and  DBHP,  statement  9  and  preced- 
ing, will  be  needed.  Statements  for  DBHP 
are  mathematical  expressions  of  the  first  part 
of  table  1,  with  diameter  as  the  dependent 
variable.  BREAK  is  the  largest  permissible 
basal  area  of  the  first  equation  and  the 
lowest  basal  area  of  the  second.  Statements 
that  contain  the  ratio  of  DLEV  or  THIN  to 


80.0  must  be  replaced  if  another  system  of 
labeling  growing  stock  levels  is  used.  This 
change  was  discussed  in  the  description  of 
table  1. 

3.  Statements  for  TOTO  and  TOTT,  so  cubic 
volumes  per  acre  are  correct  for  the  species 
and  tree  volume  equation  selected. 

4.  Statement  for  PDBHE,  based  on  program 
THIN1  and  stand  tables  for  the  species. 

5.  Statement  for  DBHO,  based  on  a  growth 
study  of  the  species  of  interest. 

6.  Table  headings. 

It  may  be  necessary  to  account  for  the  trees 
that  die  during  the  periods  between  thinnings. 
If  so,  the  statement  DENO=DENT  must  be 
replaced.  It  may  be  possible  to  relate  mortality 
in  numbers  of  trees  to  average  stand  diameter 
(Myers  1967). 


Test  Problem  for  YIELD 

A  problem  is  included  (appendix  1)  to 
provide  additional  description  of  the  data  deck 
and  of  the  output.  The  test  problem  should 
be  run  to  check  accuracy  of  punching  of  the 
source  deck  and  to  test  compatibility  with  local 
equipment. 

Values  used  in  the  test  problem  were  as 
follows: 

NTSTS  -  2  RINT  -  10.0 

JCYCL  -  20        ROTA  -  130.0 

MIX  -  5  SITE  -  60.0 

AGEO  -  30.0       THIN  -  120.0  and  80.0 

DBHO  -  4.5        TABLl(K)  -  Values  in  table  1 

DENO  -  1000.0    TABL2(K,L)  -  Values  in  table  2 

DSTY  -  70.0        TABL3(K)  -  Values  in  table  3 

PRET  -  27  0       TABL4(K)  -  Values  in  table  4 


Values  were  the  same  for  the  two  tests 
except  for  THIN.  The  test  problem  compares 
the  effect  on  yields  of  initial  thinning  to 
growing  stock  levels  80  and  120,  and  of  rethin- 
ning  to  each  of  several  levels.  One  card  or 
set  of  cards  of  each  type,  except  type  1,  was 
read  for  each  test  performed. 

Output  of  the  tests  can  be  used  in  many 
ways  to  assist  in  decisionmaking.  Yields  and 
numbers  of  noncommercial  cuts  can  be 
compared  (table  5).  Money  yields  and  rates 
earned  can  be  computed  if  necessary  data  on 
costs  and  stumpage  values  are  available.  Stand 
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Order  and  Contents  of  Data  Deck  for  Program  YIELD 


No.  of 


Card 
type 

cards 
in  type 

Var  iab 1 e 
name 

Columns 

Format 

Description  of  variable 

1 

1 

NTSTS 

•4 

14 

Number  of  tests  in  a  batch.     The  number  of  sets 
of  yield  tables  to  be  computed  and  printed  for 
each  combination  of  values  on  a  card  of  type  3. 

2 

1 

JCYCL 

1- 

■4 

14 

Interval  between  intermediate  cuts.     A  multiple 
of  RINT. 

MIX 

5- 

•8 

14 

Number  of  stocking  levels  or  values  of  DLEV  to 
be  examined  in  one  test. 

3 

1 

AGEO 

1- 

•8 

F8. 

,3 

Initial  age  to  be  shown  in  a  yield  table.  Stani 
age  when  first  thinning  occurs. 

DBHO 

9- 

•16 

F8. 

.3 

Average  stand  d.b.h.  just  prior  to  initial  thin- 
ning at  stand  age  AGEO. 

DENO 

17- 

•24 

F8. 

,3 

Number  of  trees  per  acre  just  prior  to  initial 
thinning  at  age  AGEO. 

DSTY 

25- 

•32 

F8. 

,3 

Lowest  growing  stock  level  for  intermediate  cut; 
after  initial  thinning.     Level  will  increase  by 
10  as  many  times  as  specified  by  MIX  on  card 
type  2. 

PRET 

33- 

•40 

F8. 

,3 

Estimated  percentage  of  trees  to  be  retained 
after  initial  thinning  at  age  AGEO. 

RINT 

41- 

•48 

F8. 

,3 

Number  of  years  for  which  the  growth  equation 
makes  one  projection  of  growth.     Value  is  10.0 
for  the  equation  for  DBHO. 

ROTA 

49- 

•56 

F8. 

,3 

Final  age  for  which  stand  data  are  to  be  given 
in  a  yield  table. 

SITE 

57- 

•64 

F8. 

,3 

Site  index  for  the  species. 

THIN 

65- 

•72 

F8.3 

Growing  stock  level  for  initial  thinning  at  age 
AGEO.     May  equal  DLEV. 

4 

4 

TABLl(K) 

1- 

•63 

21F3. 

1 

Table  1  of  this  publication  or  a  similar  table 
giving  basal  areas  after  thinning  in  relation 
to  average  stand  d.b.h. 

5 

3 

TABL2 (K,L) 

1- 

75 

25F3. 

1 

Table  2  of  this  publication  or  a  similar  table 

U 1.      Lice     llclgLl  Lo      U  V     Clg,  C     O.LHJ.     O  -L  L  C     111UCA  i — Laooco  i 

6 

8 

TABL3 (K) 

1- 

72 

24F3. 

3 

Table  3  of  this  publication  or  a  similar  table 
of  factors  for  conversion  of  total  cubic  feet 
to  merchantable  cubic  feet. 

7 

6 

TABL4 (K) 

1- 

78 

26F3. 

2 

Table  4  of  this  publication  or  a  similar  table 
of  factors  for  conversion  of  total  cubic  feet 
to  board  feet. 
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Table  5. — Comparison  of  alternatives  included  in  the  test  problem 


Initial  growing 

Subsequent 

Yields  per 

acre 

Numb  er 

of 

stock  level 

stocking  levels 

Intermediate 

Final 

precommercial 

thinnings 

120 


70 
80 
90 
100 
110 


Cu.  ft. 

2450 
2470 
2360 
2290 
2270 


M.  bd.  ft, 

13.8 
14.9 
15.9 
17.0 
17.7 


80 


70 
80 
90 
100 
110 


2120 
1840 
1810 
1860 
1940 


14.1 
15.2 
16.3 
17.4 
18.1 


Final  cut  in  board  feet  at  age  130.  Intermediate  cuts  in  cubic  feet.  Minimum  commercial  limit 
300  cubic  feet. 


ages  at  culmination  of  mean  annual  increment, 
and  rates  earned  assist  in  the  determination  of 
suitable  rotations.  Under  the  conditions 
specified  for  the  test  problem,  rotations  of 
about  130  years  and  initial  thinnings  that  leave 
relatively  large  growing  stocks  appear  worthy 
of  further  examination. 


Description  of  Program  THIN1 

The  FORTRAN  statement  for  PDBHE  that 
follows  statement  8  of  program  YIELD  computes 
the  logarithm  of  estimated  stand  diameter  after 
initial  thinning.  Independent  variables  are  the 
logarithms  of  diameter  before  thinning  and 
percentage  of  trees  left  after  thinning.  Equa- 
tions for  PDBHE  are  calculated  from  data  pro- 
duced by  THIN1  (appendix  2).  The  program 
simulates  several  intensities  of  thinning  from 
below  in  young  stands  with  various  average 
diameters. 

Preliminary  computations  and  operations 
performed  by  THIN1,  described  as  though  done 
manually  with  cards,  are  as  follows: 

1.  Obtain  cumulative  percentage  stand  tables 
by  1-inch  diameter  classes,  as  described  in 
item  1  of  the  section  "Information  Needed." 

2.  Convert  each  cumulative  table  to  a  stand 
table  for  1,000  trees,  and  obtain  the  average 
diameter  of  each  stand  from  average  basal 


area.  All  trees  in  a  class  are  assigned  the 
diameter  of  the  class  midpoint. 

3.  Write  the  diameter  of  each  tree  of  a  stand 
on  a  card  to  create  a  deck  of  1,000  cards  that 
represents  trees  of  a  stand. 

4.  Shuffle  the  cards  so  they  become  arranged  in 
random  order. 

5.  Select  a  percentage  of  trees  to  be  retained. 
The  1,000  cards  will  be  separated  into  groups 
of  equal  size.  One  tree  of  each  group  will 
be  left  standing  after  thinning.  Percentage 
retained  equals  1/group  size.  Vary  group 
size  (no  fractional  cards)  until  the  desired 
percentage  is  approximated  as  closely  as 
possible.  For  example,  separation  of  1,000 
cards  into  groups  of  20  cards  each  gives  a 
percentage  retention  of  1/20  or  5  percent. 

6.  Divide  the  deck  of  1,000  randomly  arranged 
cards  into  groups  of  the  desired  size,  and 
tally  the  largest  diameter  in  each  group. 
This  produces  a  postthinning  stand  table. 

7.  Compute  average  diameter  from  average 
basal  area.  Record  initial  diameter,  post- 
thinning  diameter,  and  percentage  of  trees 
retained  as  one  set  of  data  from  a  simulated 
plot.  Comparison  of  results  of  this  procedure 
with  values  from  thinned  permanent  plots 
gave  differences  of  0.1  inch  or  less. 

The  accompanying  tabulation  of  order  and 
contents  of  the  data  deck  for  THIN1  lists  the 
variables  for  which  numerical  values  are  needed. 


Order  and  Contents  of  the  Data  Deck  for  Program  THIN1 


No.  of 

Card        cards  Variable 

type      in  type  name  Columns  Format  Description  of  variable 


1  1  TERM  1-4  F4.0 

NDMS  5-8  14 

NGRP  9-12  14 

2  1  GRPS(I)  1-50  10F5.0 

3  1  AVDM  1-5  F5.1 


A  whole  number  between  0  and  1023  in  the  pseu- 
dorandom number  generator.     Change  in  the  value 
of  TERM  changes  the  starting  point  of  the  se- 
quence of  1024  pseudorandom  numbers. 

Number  of  stand  tables  to  be  tested. 

Number  of  group  sizes  to  be  tested.     This  is  the 
number  of  values  entered  for  GRPS(I)  on  card 
type  2. 

Number  of  trees  in  each  group,  from  which  one 
will  be  chosen  for  the  residual  stand.  Group 
sizes  usually  range  from  2  to  20. 

Average  d.b.h.  before  thinning  of  one  stand  to 
be  examined. 


NCLS  6-7  12  Number  of  1-inch  d.b.h.   classes  in  the  stand 

that  has  average  d.b.h.  AVDM. 

4  1  DBHC(I)  1-50  10F5.1  List  of  the  class  midpoints  of  the  1-inch  d.b.h. 

classes  in  the  stand  described  by  card  type  3. 
List  in  order  of  increasing  size. 

5  1  IFRQ(I)  1-40  1014  Cumulative  numbers  of  trees  for  the  d.b.h. 

classes  listed  on  card  type  4.     Order  of  list- 
ing corresponds  to  DBHC(I). 


Additional  description  of  the  program  is  provided 
by  the  test  problem  in  the  following  section 
and  in  appendix  2. 

Arrays  of  diameters  are  arranged  in  random 
order  in  THIN1  through  use  of  a  pseudorandom 
number  generator.  An  array  of  1,000  numbers 
is  produced  by  a  generator  of  the  form  (Green- 
bergerl961): 

Xi+1  =  AXi  +  C    (mod  2P) 

Terms  A,  C,  and  2P  are  given  constant  values 
in  statement  7  of  the  source  program.  An 
initial  value  of  X  •  is  read  in  as  variable  TERM 


so  the  series  can  be  entered  at  numerous 
points. 

Sets  of  data  produced  by  THIN1,  transformed 
as  necessary,  are  used  in  multiple  regression 
analysis  to  obtain  an  equation  for  PDBHE. 
It  may  be  necessary  to  include  more  combina- 
tions of  independent  variables  than  the  20 
shown  in  the  test  problem,  and  to  extend  their 
range  of  values.  Enough  combinations  should 
be  sampled  to  include  desired  values  of  the 
independent  variables  and  to  permit  verification 
of  transformations  of  the  variables  to  linear 
relationships.  Library  programs  for  multiple 
regression  are  available  at  most  computing 
centers. 
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Test  Problem  for  THIN1 

Values  read  in  for  the  test  problem  of 
THIN1  (appendix  2)  were  as  follows.  Note 
that  card  types  3  to  5,  inclusive,  were  included 
once  for  each  stand  table  examined. 


Decimal  points  identify  numbers  in  F-format. 
Lengths  of  fields  and  formats  of  the  numbers 
are  given  in  the  list  of  contents  of  the  data 
deck.  Diameter  class  midpoints  are  at  the 
half-inch  values  because  of  the  frequent  use  of 
full-inch  classes  for  ponderosa  pine. 


Card 

type 


Values 


1  21.0  4  5 

2  2.0  3.0  4.0  5.0  10.0 

3  2.5  5 

4  0.5  1.5  2.5  3.5  4.5 

5  165  433  749  983  1000 

3  3.0  5 

4  0.5  1.5  2.5  3.5  4.5 

5  86  265  563  877  1000 

3  3.4  6 

4  0.5  1.5  2.5  3.5  4.5  5.5 

5  25  130  408  760  969  1000 

3  4.0  6 

4  1.5  2.5  3.5  4.5  5.5  6.5 

5  11  165  574  848  966  1000 
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APPENDIX  1 
Program  YIELD  and  Output  of  Test  Problem 
Program  YIELD 


DEFINITIONS   OF   VARIABLES    IN   PROGRAM  YIELO 

AGEO  =   INITIAL   AGE    IN   YIELD  TABLE 

BASC   =   BASAL   AREA  CUT   PER  ACRE 

BASO  =   BASAL   AREA  PER   ACRE   BEFORE  THINNING 

BAST  =   BASAL   AREA  PER   ACRE   AFTER  THINNING 

BDFC  =   BOARD   FEET  CUT   PER  ACRE 

BDFO  =   BOARD  FEET  PER   ACRE   BEFORE  THINNING 

BDFT  =   BOARD  FEET  PER   ACRE   AFTER  THINNING 

CFMC  =   MERCHANTABLE   CU.FT.   CUT   PER  ACRE 

CFMO  =   MERCHANTABLE   CU.FT.   PER   ACRE   BEFORE  THINNING 

CFMT  =   MERCHANTABLE   CU.FT.    PER   ACRE   AFTER  THINNING 

DBHO  =   AVERAGE   STAND   D.B.H.    BEFORE  THINNING 

DBHT  =   AVERAGE   STAND   D.B.H.    AFTER  THINNING 

DENC   =   TREES  CUT   PER  ACRE 

DENO  =   TREES   PER   ACRE   BEFORE  THINNING 

DENT  =   TREES   PER   ACRE   AFTER  THINNING 

OLE V  =   GROWING   STOCK   LEVEL   FOR    INTERMEDIATE  CUTS 

DSTY  =   LOWEST   VALUE   OF   DLEV  USED    IN   A  TEST 

HTSO  =   TREE   HEIGHT   BEFORE  THINNING 

HTST  =  TREE   HEIGHT   AFTER  THINNING 

JCYCL  =    INTERVAL   BETWEEN    INTERMEDIATE  CUTS 

MIX  =  NUMBER  OF   STOCKING   LEVELS   EXAMINED   PER  TFST 

NTSTS   =   NUMBER   OF   TESTS   PER   BATCH, EACH  WITH   SEPARATE   HEADER  CARDS 
PRET  =   PERCENTAGE   OF   TREES   RETAINED  AFTER    INITIAL  THINNING 
R I  NT  =   NUMBER   OF   YEARS   FOR  WHICH   GROWTH  PROJECTION    IS  MADE 
ROTA  =   FINAL   AGE    IN  YIELD  TABLE 
SITE   =   SITE  INDEX 

TABL1   =   TABLE    1   OF   THIS   PUBLICATION   OR   SIMILAR  DATA 
TABL2   =   TABLE   2   OF    THIS   PUBLICATION   OR   SIMILAR  DATA 
TABL3   =   TABLE    3   OF   THIS   PUBLICATION   OR   SIMILAR  DATA 
TABL4  =   TABLE   4  OF   THIS   PUBLICATION   OR   SIMILAR  DATA 
THIN  =  GROWING   STOCK   LEVEL   FOR    INITIAL  THINNING 
TOTC   =   TOTAL   CUBIC   FEET  CUT   PER  ACRE 
TOTO  =   TOTAL   CUBIC   FEET   PER   ACRE   BEFORE  THINNING 
TOTT  =   TOTAL   CUBIC   FEET   PER   ACRE   AFTER  THINNING 

DIMENSION  TABL1 (81) ,TABL2( 4, 18 ) ,TABL3(192),TA6L4(156I 

READ  NUMBER   OF   TESTS   PER  BATCH 

READ    (5,1)  NTSTS 

1  FORMAT  (14) 

DO   31    I    =    1, NTSTS 

ZERO  VARIABLES   AND   READ    IN   VALUES   FOR    A  TEST 

BDFO  =  0.0 
BDFT  =  0.0 
CFMO  =  0.0 
CFMT  =  0.0 
DLEV   =  0.0 

READ    (5,2)    JCYCL, MIX 

2  FORMAT  (214) 

READ    (5,3)    AGEO, DBHO, DENO, DSTY, PRET, RINT, ROTA, SITE, THIN 

3  FORMAT  (9F8.3) 

READ    (5,4)    (TABLl(K) , K= 1 , 8 1 ) 

4  FORMAT  (21F3.1) 

READ    (5,5)    ( ( TABL2 ( K , L ) ,K=1 ,4) , L=l , 18 ) 

5  FORMAT  (25F3.1) 

READ    (5,6)  <TABL3(K),K=1,192) 

6  FORMAT  (24F3.3) 

READ    (5,7)    (TABL4IK) ,K=1 , 156) 

7  FORMAT  126F3.2) 
AGE01  =  AGEO 
DBH01   =  DBHO 
DEN01   =  DENO 

PROVIDE   FOR   SEVERAL   GROWING   STOCK   LEVELS   PER  TEST 
DO   31    M  =    1    ,  MIX 


OBTAIN  HTSO  FROM   TABLE  2 
TEN   =  10.0 

ISITE   =    (SITE/TEN-3.0)    +  0.01 
I  AGE □  =    (AGEO/TEN)    +  0.01 
HTSO  =   TABL2I ISITE, IAGEO) 

COMPUTE   TOTAL   CUBIC   FEET   PER  ACRE 

TOTO   =    ((0.4047   »   BASO   »   HTSO)    +    (25.5970   »   DBHO))  -191.6433 

COMPUTE   MERCH.   CU.FT.    IF   D.B.H.    IS   AT   LEAST   5.0  INCHES 


BASE)  -0.17858 


C   COMPUTE   D.B.H.    AFTER   INITIAL  THINNING 

C 

8  DO    13   J   =    1    ,  100 

PDRHE   =    (0.95462    *   ALOG 1 0 ( DBHO ) ) - ( 

DBHE   =    10.0   «»  PDBHE 

IDBHE   =   OBHE   •    10.0   +  0.5 

DBHE   =  IDBHE 

DBHE   =  DBHE/TEN 

DENE   =   DFNO   *    (PP.  ET/100.0) 

BASE   =    (0.0054542   *   DBHE   »  DBHE)  « 

BRFAK   =   49.9   *    t THIN/80.0) 

IF1BASE.GT. BREAK)    GO   TO  9 

DBHP   =    ( CO. 0/THIN)»( 0.08733   *  BASE 

GO   TO  10 

9  DBHP   =    (80.0/THIN)    «  (0.10938 
10      IOBHP  =   DBHP   »    10.0   +  0.5 

DBHP   =    I DBHP 
DBHP   =  DBHP/TEN 
I F ( DBHP-DBHE  )  11,14,12 

11  PRET   =   PRET   +   PRET   •  0.02 
GO   TO  13 

12  PRET   =   PRET   -   PRET   *  0.02 

13  CONTINUE 

14  DBHT   =  D3HE 


C   COMPUTE   BASAL   AREA   AFTER    INITIAL  THINNING 
C 

JDBHT   =    ((DBHT   -   2.0)    *    10.0)    +  1.01 

SQFT   =   TABL 1 ( JDBHT ) 

BAST   =    (THIN/80.0)    «  SQFT 

C 

C   ENTER   LOOP   FOR   REMAINING   COMPUTATIONS   AND  PRINTOUT 

C 

DO   29   K   =    1    ,  100 

IF    ( AGEO. GE. ROTA)      GO   TO  15 

DENT   =   BAST/IO. 0054542   *   DBHT   *  DBHT) 

HTST   =  HTSO 

C 

C   COMPUTE   TOTAL   CUBIC   FEET  AFTER  THINNING 
C 

TOTT   =    ((0.4047   »   BAST   *   HTST)    +    (25.5970   *   DBHT))  - 

C 

C  COMPUTE  MERCH.  CU.  FT.  IF  D.B.H.  IS  AT  LEAST  5.0  INCHES 
C 

IF    (DBHT. LT. 5.0)    GO   TO  15 

I  DBHT   =    UOBHT   -   5.0)    *    10.0)    +  1.01 

XDBHT   =   TABL3    ( IOBHT) 

CFMT   =   TOTT    *  XDBHT 


ALOGIO(PRET) 1+0.26959 


COMPUTE   BD.    FT.    IF  D. 


IS   AT  LEAST 


IF    (D6HT.LT.8.0)    GO   TO  15 

JDBHT   =    ((DBHT   -   8.0)    »    10.0)    +  l.C 

YDBHT   =   TABL4 ( JDBHT ) 

BDFT   =   TOTT   *  YDBHT 

C 

C  CHANGE  MODE  AND  ROUND  OFF  FOR  PRINTING 
C 

15  JAGEO  =  AGEO 
JSITE  =  SITE 
JDENO  =  DENO 
JHTSO 
JTOTO 
JTOTO 
JBASO 
JCFMO 
JCFMO 
JBDFO 
JBDFO 
JDENT 
JHTST 
JTOTT 
JTOTT 
JCFMT 
JCFMT 


HTSO 

(TOTO/TEN)    +  0.5 
JTOTO  *  10 
BASO  +  0.5 
(CFMO/TEN)    +  0.5 
JCFMO  ♦  10 
(BDFO/100.0)    +  0.5 
JBDFO   «  100 
DENT 
HTST 

( TOTT/ T  EN )    +  0.5 
JTOTT   ♦  10 
(CFMT/TEN)    +  0.5 
JCFMT   *  10 
IF    I JCFMT. GT. JCFMO)    JCFMO   =  JCFMT 
JBDFT   =    (BDFT/100.0)    +  0.5 
JBDFT   =   JBDFT   *  100 

IF    ( JBDFT. GT. JBDFO)    JBDFO  =  JBDFT 

JBAST   =   BAST  +  0.5 

JDENC   =   JDENO  -  JDENT 

JBASC   =    JBASO  -  JBAST 

JTOTC   =   JTOTO  -  JTOTT 

JCFMC   =   JCFMO  -  JCFMT 

JBDFC   =   JBDFO  -  JBDFT 


IF    (DBHO. LT. 5.0)    GO  TO  8 
I DBHO  =    ( (DBHO-5.0)»10.01 
XDBHO   =   TABL3I I DBHO 1 
CFMO  =   TOTO   »  XDBHO 


1.01 


COMPUTE   BD.    FT.    IF   D.B.H.    IS   AT   LEAST   8.0  INCHES 
TO  8 

01 


C   WRITE   HEADINGS   FOR   YIELD  TABLE 

C 

IF    (K.GE.2)    GO   TO  22 
WRITE    (6,17)    JSITE, JCYCL 

17  FORMAT  (1H1/1H-,27X,81HYIELDS  PER  ACRE  OF  MANAGED,  EVEN-AGED  STAND 
IS  OF  PONDEROS A  PINE  IN  THE  BLACK  HILLS/1H  ,49X,10HSITE  INDEX, 13, 1H 
2,,    14,    19H-YEAR  CUTTING  CYCLE) 

WRITE  (6,18) 

18  FORMAT    ( 1H0,25X,38HENTIRE    STAND  BEFORE   AND   AFTER   THINNING, 28X,26HP 
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IERIODIC   CUT   AND  MORTALITY) 
WRITE  (6,19) 

19  FORMAT  (1  HO, 9X.5HSTA NO, 10X,5H8ASAL,3X,7HAVERAGE,2X,7H AVERAGE, 3X.5H 
1T0TAL,3X,9HMERCHANT-,3X,9HSAWTIMBER,9X,5H8ASAL,4X,5HT0TAL,3X,9HMGR 
2CHANT-.3X.9HSAWTIMBER) 

WRITE  (6,20) 

20  FORMAT  (1H  ,  1  OX  ,  3H  AGE  ,  4X  ,  5HTRE  E  S  ,  3X  ,  4H  AR  E  A  ,  4X  ,  6HD  .  0  .  H  .  ,  3X  ,  6HHE  I  GHT 
1,2X,6HV0LUME,2X,11HABLE  VOLUME , 4X , 6H VOLUME , 3X , 5HTREES , 3X , 4H AREA , 3X 
2,6HV0LUKE,?X,11H,\BLE   VOLUME, 4X,6H VOLUME ( 

WRITE  (6,21) 

21  FORMAT  (1H  , 8X, 7H ( YEARS) , 3X.3HNO. ,3X,6HSQ. FT. , 4X,3HIN. , 6X, 3HFT. , 4X 
1 , 6HCU.FT. , 5X .6HCU.FT. ,6X,6HBD.FT.,4X,3HN0.,3X,6HSQ.FT.,2X,6HCU.FT. 
?,5X,6HCU.FT.,6X,6HBD.FT. ) 

22  WRITE    (6,23)    J AGEO , J DENO , J BA SO , ORHO , JHTSO , JTOTO , JC FMO , J BDFO 

2  3   FORMAT  (1H0,9X,I4,4X,I5,2X,I4,5X,F5.1,5X,I3,4X,I5,6X,I5,6X,I6) 
IF    (AGEO. GE. ROTA)   GO  TO  30 

WRITE    (6,24)    J AGEO, JOENT, JBAST.DBHT, JHTST, JTOTT, JCFMT , JBDFT, JDENC, 
1JBASC, JTOTC.JCFMC, JBDFC 
24   FORMAT    (1H    ,9X, I4.4X, I5.2X, I4.5X, F5. 1,5X, I 3,4X, 15, 6X, I5.6X, I6,4X , [ 
15,3X,I3,5X,I4,6X,I4,8X,I5) 

C 

C   COMPUTE   VALUES   FOR   EACH   PERIOD.    THIN   AS  SPECIFIED 
C 

I R I  NT   =   R I  NT 

IK   =  JCYCL/IRINT 

DO   26   L   =  1,IK 

AGFO  =  AGEO  +  R I  NT 

IF    ( AGEO. GT. ROTA)    GO   TO  30 

c 

C   COMPUTE   NEW   D.B.H.    BEFORE   THINNING   AND  ROUND   OFF   TO   0.1  INCH 
C 

DBHO  =    1.0097   »   DBHT   +   0.0096   *   S I T E- ( 1 . 5766* ALOG 1 0 ( 8 AS T ) ) + 3 . 302 1 
IDBHO   =   DBHO   *    10.0  +  0.5 
DBHO  =  IDBHO 
DBHO  =  OBHO/TEN 
DENO  =  DENT 

BASO  =   DENO   *(0.0054542    *   DBHO   ♦  DBHO) 

C 

C   OBTAIN  HTSO   FROM   TABLE  2 
C 

ISITE   -    ( SITE/TEN-3.0)    +-  0.01 
I  AGE  0   =    (AGEO/TEN)    +  0.01 
HTSO  =  TABL2(ISITE,IAGE0) 

C 

C  COMPUTE   TOTAL   CUBIC  FEET  BEFORE  THINNING 
C 

TOTO   =    (0.4047   *   BASO   *   HTSO)    +    (25.5970   *   DBHO ) - 1 9 1 . 643 3 

C 

C   COMPUTE   MERCH   CU.FT.    IF   D.B.H.    IS   AT   LEAST    5.0  INCHES 
C 

IF    (DBHO. LT. 5.0)    GO   TO  25 


IDBHO   =    ( (DBH0-5.0)»10.0)+1.01 
XD8H0   =   T ABL3 ( IDBHO) 
CFMO  =   TOTO   *  XDBHO 

C 

C   COMPUTE   BD.FT.IF   DBH   IS   AT  LEAST   8.0  INCHES 

C 

IF    (D8H0.LT. 8.0)    GO   TO  25 
JD8HO  =    ( (DBHO-8. 01*10.0)    +  1.01 
YDBHO  =   T A8L4 ( J  DBHO ) 
BDFO  =   TOTO   *  YDBHO 

25  IF    (L.EQ.IK)    GO   TO  27 

C 

C      WRITE   VALUES   FOR   PERIOD.    IF   THINNING  NOT  DUE 

C 

KDENO  =  DENO 

KAGEO  =  AGEO 

KHTSO  =  HTSO 

KB A  SO  =   BASO   +  0.5 

KTOTO  =    (TOTO/TEN)    +  0.5 

KTOTO  =   KTOTO   *  10 

KCFMO  =    (CFMO/TEN)    +  0.5 

KCFMO  =   KCFMO   *  10 

KBDFO  =    (BDFO/100.0)    +  0.5 

KBDFO  =   KBDFO   *  100 

WRITE    (6,23)    KAGEO, KDENO, KBASO, DBHO, KHTSO, KTOTO, KCFMO, KBDFO 
DBHT   =  DBHO 
BAST   =  BASO 

26  CONTINUE 

C 

C  INCREASE   D.B.H.    BY   THINNING   AND  COMPUTE   POST-THINNING  VALUES 

C 

27  DBHT   =   DBHO   +  0.4 

IF    (DBHT.GE.1C.0)    GO   TO  28 

J  DBHT   =    (  (DBHT-2.0)*10.0)  +1.01 

SQFT   =   TA6L 1 ( JDBHT ) 

BAST   =  (DLEV/80.0)*SQFT 

GO  TO  29 

28  3AST   =  DLEV 

29  CONTINUE 

C 

C   PREPARE   FOR  NEXT   TABLE   OF   THE  TEST 

C 

30  AGEO  =  AGEO  1 
BDFO  =  0.0 
SOFT  =  0.0 
CFMO  =  0.0 
CFMT  =  0.0 
DBHO  =  DBHO 1 
DENO  =  DENO 1 

31  CONTINUE 
CALL  EXIT 
END 


YIELD  Output 


YIELDS    PER   ACRE   OF   MANAGED,    EVEN-AGED   STANDS   OF    PONDEROSA   PINE    IN   THE    BLACK  HILLS 
SITE    INDEX   60,      20-YEAR  CUTTING  CYCLE 


ENTIRE    STAND   BEFORE   AND   AFTER  THINNING 


PERIODIC   CUT    AND  MORTALITY 


STAND 
AGE 
( YEARS) 

TREES 
NO. 

BASAL 
AREA 
SQ.FT. 

AVERAGE 
D.B.H. 
IN. 

AVERAGE 
HEIGHT 
FT. 

TOTAL 
VOLUME 
CU.FT. 

MERCHANT- 
ABLE VOLUME 
CU.FT. 

SAWTIMBER 
VOLUME 
BD.FT. 

TREES 
NO. 

BASAL 
AREA 
SQ.FT. 

TOTAL 
VOLUME 
CU.FT. 

MERCHANT- 
ABLE VOLUME 
CU.FT. 

SAWTIMBER 
VOLUME 
BD.FT. 

30 
30 

1000 
496 

110 
73 

4.5 

5.2 

?0 
20 

820 
530 

200 
200 

0 
0 

504 

37 

290 

0 

0 

40 

496 

104 

6.2 

28 

1150 

680 

0 

50 

50 

496 
200 

133 
60 

7.0 
7.4 

35 
35 

1870 
850 

1380 
670 

0 
0 

296 

73 

1020 

710 

0 

60 

200 

79 

8.5 

41 

1340 

1170 

1500 

70 

70 

200 
130 

99 
70 

9.5 

9.9 

47 
47 

1930 
1390 

1750 
1280 

3500 
2900 

70 

29 

540 

470 

600 

80 

130 

86 

11.0 

52 

1910 

1780 

5500 

90 
90 

130 
84 

101 
70 

11.9 
12.3 

57 
57 

2440 
1740 

2300 
1640 

8500 
6400 

46 

31 

700 

660 

2100 

100 

83 

13.4 

60 

2170 

2050 

8800 

110 
110 

84 
58 

96 

70 

14.4 
14.8 

63 
63 

2620 
1970 

2490 
1880 

11500 
8800 

26 

26 

650 

610 

2700 

120 

58 

81 

15.9 

66 

2370 

2270 

11200 

130 

58 

91 

16.9 

69 

2790 

2670 

13800 
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YIELDS   PER   ACRE   OF   MANAGED,    EVEN-AGED   STANDS   OF   PONDEROSA   PINE    IN   THE   B L AC 
SITE    INDEX   60,      20-YEAR   CUTTING  CYCLE 


ENTIRE   STAND   BEFORE    AND   AFTER  THINNING 


CUT   AND  MORTALITY 


STAND 
AGE 
(YEARS) 


TREES 
NO. 


BASAL 
AREA 
SO. FT. 


AVERAGE 
D.B.H. 
I.N. 


AVERAGE 
HE  I GHT 


TOTAL 
VOLUME 
CU.FT. 


MERCHANT- 
ABLE VOLUME 
CU.FT. 


SAWTIMBER 
VOLUME 


TREES 
NO. 


BASAL 
AREA 
SQ.FT. 


TOTAL  MERCHANT- 
VOLUME  A8LF  VOLUME 
CU.FT.  CU.FT. 


SAWTIMBER 
VOLUME 
BD.FT. 


1000 
496 


4.5 

5.2 


820 
530 


200 
200 


496 
229 


229 
152 


90 
110 


7.0 
7.4 


1870 
970 


2150 
1580 


1380 
760 


1950 
1450 


3800 
3200 


90 
100 


11.7 
12.  1 


2730 
1960 


2560 
1850 


9100 
7100 


110 
110 


120 
130 


100 
70 


14.0 
14.4 

15.4 

16.3 


2900 
2220 


2650 
3090 


2750 
2110 


2530 
2960 


12300 
9700 


12200 
14900 


I  ELDS   PER   ACRE   OF   MANAGED,    EVEN-AGED   STANDS   OF   PONDEROSA   PINE    IN   THE  BLACK 
SITE    INDEX   60,      20-YEAR  CUTTING  CYCLE 


ENTIRE   STAND   BEFORE    AND   AFTER  THINNING 


PERIODIC   CUT    AND  MORTALITY 


STAND 
AGE 
(YEARS) 


BASAL 
AREA 
SQ.FT. 


AVtRAGE 
D.B.H. 
IN. 


TOTAL 
VOLUME 
CU.FT. 


MERCHANT- 
ABLE VOLUME 
CU.FT. 


SAWTIMBER 
VOLUME 
BD.FT. 


BASAL 
AREA 
SO. FT. 


TOTAL  MERCHANT- 
VOLUME  ABLE  VOLUME 
CU.FT.  CJ.FT. 


SAWTIMBER 
VOLUME 
BD.FT. 


1000 
496 


4.5 
5.2 


820 
530 


200 
200 


496 
257 


104 
133 


1870 
1090 


1380 
860 


257 
177 


99 
119 


9.2 
9.6 


2310 
1750 


2070 
1600 


90 
100 


110 
110 


120 
120 
120 


123 
90 


11.7 
12.6 


13.4 


2950 
2180 


3160 
2460 


2760 
2050 


2990 
2330 


9000 
7300 


12900 
10300 


120 
130 


102 
115 


14.7 
15.6 


2910 
3420 


2770 
3270 


12900 
15900 


YIELDS   PER   ACRE   OF   MANAGED,    EVEN-AGED   STANDS   OF   PONDEROSA   PINE    IN   THE   BLACK  HILLS 


SITE   INDEX  60, 
ENTIRE   STAND   BEFORE   AND  AFTER  THINNING 


20-YEAR   CUTTING  CYCLE 


PERIODIC   CUT   AND  MORTALITY 


STAND 
AGE 
(YEARS) 


1000 
496 


286 
200 


BASAL 
AREA 
SQ.FT. 

110 
73 

104 

133 
86 

108 

129 


AVERAGE  AVERAGE 


D.  B. 
IN. 


4.5 
5.2 


7.0 
7.4 


9.  1 
9.5 


HEIGHT 
FT. 

20 
20 

28 

35 
35 

41 


TOTAL 
VOLUME 
CU.FT. 

820 
530 


1870 
1210 


2500 
1930 


MERCHANT- 
ABLE VOLUME 
CU.FT. 

200 
200 


1380 
950 


SAWTIMBER 
VOLUME 
BD.FT. 


TREES 
NO. 


BASAL 
AREA 
SQ.FT. 


TOTAL 
VOLUME 
CU.FT. 


MERCHANT- 
ABLE VOLUME 
CU.FT. 


SAWTIMBER 
VOLUME 
BD.FT. 


2240 
1760 


3900 
3500 


200 
141 


132 
100 


10.3 
11.0 


3140 
2410 


2940 
2260 


110 
110 


141 
102 


130 
100 


13.0 
13.4 


3460 
2700 


3260 
2560 


120 
130 


102 
102 


114 
127 


14.3 
15.  1 


3220 
3740 


3060 
3560 


14000 
17000 
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YIELDS   PER   ACRE   OF    MANAGED,    EVEN-AGED   STANDS   OF    PONDEROS A   PINE    IN   THE   BLACK  HILLS 
SITE    INDEX   60,      20-YEAR   CUTTING  CYCLE 


ENTIRE   STAND   BEFORE    AND  AFTER  THINNING 


PERIODIC   CUT   AND  MORTALITY 


STAND 

BASAL 

AVERAGE 

AVERAGE 

TOTAL 

MERCHANT- 

SAWTIMBER 

RASAL 

TOTAL 

MERCHANT- 

AGE 

TREES 

AREA 

HEIGHT 

VOLUME 

ABLE  VOLUME 

VOLUME 

TREES 

AREA 

VOLUME 

ABLE  VOLUME 

( YEARS ) 

NO. 

SQ.FT. 

IN. 

FT  . 

CU . FT . 

CU.FT. 

8D.FT. 

NO. 

SQ.FT. 

CU.FT. 

CU.FT. 

30 

1000 

1  10 

4.  5 

20 

820 

200 

0 

30 

496 

73 

5.2 

20 

530 

200 

0 

504 

37 

290 

0 

40 

496 

104 

6.2 

28 

1150 

680 

0 

50 

496 

133 

7.0 

35 

1870 

1380 

0 

50 

314 

94 

7.4 

35 

1330 

1050 

0 

182 

39 

540 

330 

60 

314 

115 

8.  2 

41 

1930 

1660 

1800 

70 

314 

136 

8.9 

47 

2620 

2330 

3700 

70 

228 

108 

9.3 

47 

2100 

1890 

3500 

86 

28 

520 

440 

80 

228 

127 

10.1 

52 

2750 

2550 

6200 

90 

228 

146 

10.8 

57 

3440 

3220 

9400 

90 

160 

no 

11.2 

57 

2630 

2460 

7900 

68 

36 

810 

760 

100 

160 

126 

12.0 

60 

3180 

2990 

11300 

110 

160 

141 

12.7 

63 

3740 

3530 

14400 

110 

117 

110 

13.  1 

63 

2950 

2790 

11700 

43 

31 

790 

740 

120 

117 

124 

13.9 

66 

3470 

3300 

14700 

130 

117 

137 

14.6 

69 

4000 

3800 

17700 

SAWTIMBER 
VOLUME 
8D.FT. 


YIELDS   PER   ACRE   OF   MANAGED,    EVEN-AGED   STANDS   OF    PONDEROSA   PINE    IN   THE   BLACK  HILLS 
SITE   INDEX  60,      20-YEAR  CUTTING  CYCLE 


ENTIRE    STAND   BEFORE   AND   AFTER  THINNING 


PERIODIC   CUT   AND  MORTALITY 


STAND 
AGE 
(YEARS) 

TREES 
NO. 

BASAL 
AREA 
SQ.FT. 

AVERAGE 
D.B.H. 
IN. 

AVERAGE 
HEIGHT 
FT. 

TOTAL 
VOLUME 
CU.FT. 

MERCHANT- 
ABLE VOLUME 
CU.FT. 

SAWT IMBER 
VOLUME 
BD.FT. 

TREES 
NO. 

BASAL 
AREA 
SQ.FT. 

TOTAL 
VOLUME 
CU.FT. 

MERCHANT- 
ABLE VOLUME 
CU.FT. 

30 
30 

1000 
320 

110 

51 

4.5 
5.4 

20 
20 

820 
360 

150 
150 

0 
0 

680 

59 

460 

0 

40 

320 

76 

6.6 

28 

840 

570 

0 

50 
50 

320 
181 

101 
63 

7.6 

8.0 

35 
35 

1430 
910 

1160 
770 

700 
700 

139 

38 

520 

390 

60 

181 

82 

9.  1 

41 

1400 

1260 

2200 

70 
70 

181 

99 
70 

10.0 
10.4 

47 
47 

1950 
1410 

1810 
1310 

4200 
3400 

63 

29 

540 

500 

80 

118 

86 

11.5 

52 

1900 

1790 

6100 

90 

90 

78 

100 
70 

12.4 
12.8 

57 
57 

2420 
1750 

2280 
1650 

9000 
6800 

40 

30 

670 

630 

100 

78 

83 

13.9 

60 

2170 

2060 

9200 

110 
110 

78 
54 

95 
70 

14.9 
15.3 

63 
63 

2610 
1980 

2490 
1890 

11700 
9100 

24 

25 

630 

600 

120 

54 

80 

16.4 

66 

2380 

2270 

11500 

130 

54 

91 

17.4 

69 

2780 

2670 

14100 

SAWT IMBER 
VOLUME 


YIELDS   PER   ACRE   OF   MANAGED,    EVEN-AGED   STANDS   OF    PONDEROSA   PINE    IN    THE    BLACK  HILLS 
SITE    INDEX  60,      20-YEAR  CUTTING  CYCLE 


ENTIRE   STAND   BEFORE   AND   AFTER  THINNING 


PERIODIC   CUT    AND  MORTALITY 


STAND 
AGE 
(YEARS) 

TREES 
NO. 

BASAL 
AREA 
SQ.FT. 

AVERAGE 
D.B.H. 
IN. 

AVERAGE 
HEIGHT 
FT. 

TOTAL 
VOLUME 
CU.FT. 

MERCHANT- 
ABLE VOLUME 
CU.FT. 

SAWT I MBER 
VOLUME 
BD.FT. 

TREES 
NO. 

BASAL 
AREA 
SQ.FT. 

TOTAL 
VOLUME 
CU.FT. 

MERCHANT- 
ABLE VOLUME 
CU.FT. 

30 
30 

1000 
320 

110 
51 

4.5 

5.4 

20 
20 

820 
360 

150 
150 

0 
0 

680 

59 

460 

0 

40 

320 

76 

6.6 

28 

840 

570 

0 

50 
50 

320 
207 

101 
73 

7.6 
8.0 

35 
35 

1430 
1040 

1160 
880 

800 
800 

113 

28 

390 

280 

60 

207 

92 

9.0 

41 

1560 

1390 

2300 

70 
70 

207 
136 

111 

80 

9.9 
10.  3 

47 
47 

2170 
1590 

2010 
1480 

4600 
3800 

69 

31 

580 

530 

80 

138 

96 

11.3 

52 

2120 

1990 

6500 

90 
90 

138 

92 

112 

80 

12.2 
12.6 

57 
57 

2710 
1980 

2550 
1860 

9900 
7500 

46 

32 

730 

690 

100 

92 

93 

13.6 

60 

2420 

2290 

10000 

110 
110 

92 
66 

106 

80 

14.5 
14.9 

63 
63 

2880 
2230 

2740 
2120 

12600 
10000 

26 

26 

650 

620 

120 

66 

91 

15.9 

66 

2650 

2530 

12500 

130 

66 

102 

16.8 

69 

3080 

2950 

15200 

SAWT IMBER 
VOLUME 
3D. FT. 
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YIELDS   PER   ACRE   OF   MANAGED,    EVEN-AGED   STANDS   OF   PONDEROSA   PINE    IN   THE  BLACK 
SITE    INDEX   60,      20-YEAR   CUTTING  CYCLE 


ENTIRE   STAND   BEFORE    AND   AFTER  THINNING 


PERIODIC   CUT   AND   MORTAL  I TY 


STAND 
AGE 
(YEARS) 


BASAL 
AREA 
SQ.FT. 


AVERAGE 
D.B.H. 
IN. 


AVERAGE 
HE IGHT 


TOTAL 
VOLUME 
CU.FT. 


MERCHANT- 
ABLE VOLUME 
CU.FT. 


SAWT I MBEK 
VOLUME 
BD.FT. 


BASAL 
AREA 
SQ.FT. 


TOTAL 
VOLUME 
CJ.FT. 


MERCHANT- 
ABLE VOLUME 
CU.FT. 


SAWTIMBcR 
VOLUME 
BD.FT. 


320 
233 


7.6 
8.0 


1430 
1170 


570 
L160 


900 
900 


90 
100 


110 
110 


233 
161 


161 
110 


110 
110 


120 
90 


123 
90 


104 
117 


11.8 
12.2 


13.9 
14.3 


2340 
1780 


2940 
2200 


3140 
2470 


2140 
1650 


2760 
2070 


4600 
4000 


6700 
10100 


120 
130 


102 
114 


15.2 
16.1 


2910 
3410 


2780 
3260 


13300 
16300 


YIELDS   PER   ACRE   OF   MANAGED,    EVEN-AGED   STANDS   OF   PONDEROSA   PINE    IN   THE   BLACK  HILLS 
SITE    INDEX   60,      20-YEAR  CUTTING  CYCLE 


ENTIRE   STAND   BEFORE   AND   AFTER  THINNING 


CUT   AND  MORTALITY 


STAND 
AGE 
(YEARS) 


BASAL 
AREA 
SQ.FT. 


AVERAGE 
D.B.H. 
IN. 


AVERAGE 
HEIGHT 


TOTAL 
VOLUME 
CU.FT. 


MERCHANT- 
ABLE VOLUME 
CU.FT. 


SAWTIMBER 
VOLUME 
BD.FT. 


BASAL 
AREA 
SQ.FT. 


TOTAL 
VOLUME 
CU.FT. 


MERCHANT- 
ABLE VOLUME 
CU.FT. 


SAWTIMBER 
VOLUME 
BD.FT. 


110 
51 


320 
259 


6.6 
7.6 


1430 
1300 


259 
183 


130 
100 


9.6 
10.0 


2540 
1970 


2320 
1820 


4800 
4300 


183 
129 


132 
100 


10.8 
11.5 


3150 
2420 


2960 
2270 


10100 
8400 


110 
110 


129 
129 


129 
100 


3.5 
3.9 


3440 
2710 


3250 
2580 


14100 
11500 


120 
130 


113 

126 


3220 
3720 


3060 
3560 


14400 
17400 


YIELDS    PER   ACRE   OF   MANAGED,    EVEN-AGED   STANDS   OF   PONDEROSA   PINE   IN   THE   BLACK  HILLS 
SITE    INDEX   60,      20-YEAR   CUTTING  CYCLE 


ENTIRE   STAND   BEFORE    AND   AFTER  THINNING 


PERIODIC   CUT   AND  MORTALITY 


STAND 
AGE 
(YEARS) 


BASAL 
AREA 
SQ.FT. 


AVERAGE 
D.B.H. 
IN. 


AVERAGE 
HEIGHT 


TOTAL 
VOLUME 
CU.FT. 


MERCHANT- 
ABLE VOLUME 
CU.FT. 


SAWTIMBER 
VOLUME 
BD.FT. 


BASAL 
AREA 
SQ.FT. 


TOTAL 
VOLUME 
CU.FT. 


MERCHANT- 
ABLE VOLUME 
CU.FT. 


SAWTIMBER 
VOLJME 
BD.FT. 


1000 
320 


110 
51 


5.4 
6.6 


820 
360 


150 
150 


320 
285 


101 
100 


7.6 
8.0 


285 
205 


121 
141 


9.5 
9.9 


2730 
2150 


2480 
1990 


146 
110 


3460 
2650 


3240 
2490 


10900 
9100 


140 
110 


13.3 
13.7 


120 
130 


107 
107 


123 
135 


14.5 
15.2 


3470 
3980 


3300 
3800 


15200 
18100 
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APPENDIX  2 


Program  TfflNl  and  Output  of  Test  Problem 
Program  TfflNl 


c 

11 

FORMAT  (F5.1.I2) 

c 

DEFINITIONS   OF  VARIABLES 

ADBU(N)    =  AVDM 

c 

ADBO 

AVERAGE    DIAMETER   BEFORE  THINNING 

READ    (5,12)    (DBHCI I ) , 1=1 ,NCLS ) 

c 

ADBT 

AVERAGE   DIAMETER    AFTER  THINNING 

12 

FORMAT  (10F5.1) 

c 

AVBA 

- 

AVERAGE   BASAL   AREA  OF  RESERVED  TREES 

READ    (5,13)    (IF RO (I), 1  =  1, NCL  S ) 

c 

AVDM 

AVERAGE    STAND   DBH   BEFORE  THINNING 

13 

FORMAT  (1014) 

c 

BAST 

BASAL    AREA   OF   RESERVED  TREES 

C 

c 

DBHC 

= 

DIAMETER   CLASS   OF   A   STAND  TABLE 

C 

EXPAND   STAND   TABLE   TO   1000    INDIVIDUAL  DIAMETERS 

c 

D  I  AM 

= 

INDIVIDUAL   DIAMETERS   FROM   DBHC   AND  IFRQ 

C 

c 

IFRQ 

CUMULATIVE   NUMBER   OF   TREES   BY   DBH  CLASSES 

IX  =  1 

c 

IRND 

RANDOM  NUMBERS   FROM    1    TO  1000 

DO    15  K=1,NCLS 

c 

NCLS 

NUMBER   OF   ONE-INCH   DBH   CLASSES    IN   STAND  TABLE 

JX   =  IFRQ(K) 

c 

NDMS 

NUMBER   OF   AVERAGE   STAND   DIAMETERS  TESTED 

DO    14  L=IX,JX 

c 

PRE  T 

PERCENTAGE   OF   TREES   LEFT   AFTER  THINNING 

14 

DIAM(L)    =  DBHC(K) 

c 

RSRV 

DIAMETER   DF   A   RESERVED  TREE 

15 

IX   =    1    +  JX 

c 

SDSQ 

SUM  OF   SQUARED  DIAMETERS 

C 

c 

TERM 

TERM   IN   PSEUDORANDOM  NUMBER  GENERATOR 

C 

ARRANGE    INDIVIDUAL    DIAMETERS   AT  RANDOM 

c 

TREE 

INDIVIDUAL   DIAMETERS    ARRANGED   AT  RANDOM 

C 

c 

GRPS 

GROUP  SIZE   FOR  THINNING  INTENSITY 

DO   16  IK=1,1000 

c 

NGRP 

NUMBER   OF   GROUP    SIZES  TESTED 

N3R   =    I RND ( IK) 

c 

NOGP 

NUMBER   OF   GROUPS, EACH   OF   SIZE  GRPS 

DMTR   =   DIAM  (IK) 

c 

16 

TREE ( NBR I    =  DMTR 

DIMENSION  DBHC (10),IFRQ(10I,DIAM(1000),TREF(1000),IRND(1000) 
1RSRVI  1000) , PRETI 10) , ADBO I  10) , ADBT 1 10, 10)  ,GRPS(  10) 

INITIALIZE    AND   READ   VARIABLES   USED   FOR   ALL    STAND  TABLES 

DO   1  1=1,10 
ADBO ( I )    =  0.0 
GRPSI I  )   =  0.0 

1  PRETI I )    =  0.0 
DO   2  1=1,10 
DO  2  J=l,10 

2  ADBT ( I, J)   =  0.0 
00   3  1=1,1000 

3  I RND (  I  )    =  0 

READ    (5,4)    TERM, NDMS, NGRP 

4  FORMAT    (F4. 0,214) 

READ    (5,5)    ( GRPSI I ), 1=1, NGRP) 

5  FORMAT  (10F5.0) 

COMPUTE   PERCENTAGE   OF   TREES   RETAINED   FOR   EACH   GROUP  SIZE 

DO   6    1=1, NGRP 

6  PRETI  I  )    =   100.0/GRPSI I ) 


GENERATE    1000   PSEUDORANDOM  NUMBERS 


3.01/1024.0 

3)-    1024   *  NDIV 


DO   8  1=1,1000 

7  NDIV   =    117.0  *  TERM 
NTERM   =  TERM 
NTERM  =    (17   *  NTERM 
TERM  =  NTERM 

IF    (NTERM    .LE.    0)    GO   TO  7 
IF    (NTERM   .GT.    1000)    GO  TO  7 

8  IRNDI I )   =  NTERM 

DO  COMPUTATIONS   FOR  EACH  STAND  TABLE 
DO   21  N=1,NDMS 

INITIALIZE    AND   READ   FOR   EACH   STAND  TABLE 

DO   9  1=1,10 
DBHC ( I )    =  0.0 

9  IFRQ! I )   =  0 

DO   10  1=1,1000 
DIAMI I )    =  0.0 
RSRVI I  )   =  0.0 
10  TREE  I  I  )   =  0.0 

READ   (5,11)    AVDM, NCLS 


CREATE  GROUPS  OF  SPECIFIED  SIZES. 
RECORD  LARGEST   DBH   IN  EACH  GROUP 

DO  21  M=1,NGRP 
BAST  =  0.0 
AVBA   =  0.0 
SDSQ  =  0.0 

NOGP  =  1000.0/GRPSIM) 
ANOGP  =  NOGP 


IGRP  =  GRPS(M) 
MX  =  IGRP 
DO   19  IJ=1,N0GP 
DO   17  IL=1,NCLS 
JK  =  NCLS  -  IL 
DO    17  NK=LX,MX 
IF    (TREEINK)  .E 

17  CONTINUE 

18  RSRV(IJ)   =   DBHC ( JK ) 
LX   =   MX   +  1 

19  MX  = 


DBHC ( JK ) )    GO  TO 


IGRP 


COMPUTE    POST-THINNING   STAND   DIAMETER   FROM  BASAL  AREA 

DO  20  KL=1,N0GP 

20  SDSQ  =   SDSQ  +  RSRV(KL)    *  RSRVI K L ) 
BAST   =   0.0054542   •  SDSQ 

AVBA  =  BAST/ ANOGP 
DMSO  =   AVBA/O. 0054542 
ADBTIN.M)    =  SQRT(DMSQ) 

21  CONTINUE 
WRITE  (6,22) 

22  FORMAT    I 1H1 ,////, 43X,48HAVERAGE    STAND   DIAMETER   AFTER   THINNING  FROM 
1  BELOW) 

WRITE  (6,23) 

23  FORMAT    I 1H   , 54X , 26HBLACK  HILLS   PONDEROSA  PINE) 
WRITE  (6,24) 

24  FORMAT    ( 1 HO , 1 3X , 8HD I AMET ER / 1 H    , 14X , 6HBE FORE / 1H   . 1 3X , 8HTH I NN  I  NG ,  32X 
1.25HPERCENT   OF  TREES  RETAINED) 

WRITE    (6,25)    ( PRETI I ) , 1=1, NGRP) 

25  FORMAT    ( 1H    , 1 3X , 8H I  I NC HE S ) ,    6X , F4. 1 , 9 ( 7X , F4. 1 ) ) 
DO  27  MN= 1 , NDMS 

WRITE    (6,26)    ADB0(MN),(ADBT(MN,MK),MK=1, NGRP ) 

26  FORMAT    (    1H0,    16X,    F4.1,    7X,    F4.1,  9(7X,F4.1)> 

27  CONTINUE 
CALL  EXIT 
END 


THINl  Output 


AVERAGE   STAND  DIAMETER   AFTER  THINNING  FROM  BELOW 
BLACK  HILLS   PONDEROSA  PINE 

DIAMETER 
BEFORE 

THINNING  PERCENT  OF   TREES  RETAINED 


(INCHES)  50.0  33.3  25.0  20.0  10.0 

2-  5  2.9  3.1  3.3  3.4  3.6 
3.0  3.4  3.7  3.9  4.0  4.3 

3-  4  3.9  4.2  4.3  4.4  4.7 
4.0  4.6  4.9  5.1  5.3  5.7 


—  csu.  Ft.  conb. 
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